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Abstract 

We study the non-perturbative behavior of two versions of the QCD effective charge, one ob- 
tained from the pinch technique gluon self-energy, and one from the ghost-gluon vertex. Despite 
their distinct theoretical origin, due to a fundamental identity relating various of the ingredients 
appearing in their respective definitions, the two effective charges are almost identical in the en- 
tire range of physical momenta, and coincide exactly in the deep infrared, where they freeze at 
a common finite value. Specifically, the dressing function of the ghost propagator is related to 
the two form factors in the Lorentz decomposition of a certain Green's function, appearing in a 
variety of field-theoretic contexts. The central identity, which is valid only in the Landau gauge, is 
derived from the Schwinger-Dyson equations governing the dynamics of the aforementioned quan- 
tities. The renormalization procedure that preserves the validity of the identity is carried out, and 
various relevant kinematic limits and physically motivated approximations are studied in detail. A 
crucial ingredient in this analysis is the infrared finiteness of the gluon propagator, which is inextri- 
cably connected with the aforementioned freezing of the effective charges. Some important issues 
related to the consistent definition of the effective charge in the presence of such a gluon propagator 
are resolved. We finally present a detailed numerical study of a special set of Schwinger-Dyson 
equations, whose solutions determine the non-perturbative dynamics of the quantities composing 
the two effective charges. 

PACS numbers: 11.15.Tk 12.38.Lg, 12.38.Aw, 



I. INTRODUCTION 



The infrared behavior of the QCD effective charge is of considerable theoretical and 
phenomenological interest This quantity, when correctly defined, provides a 

continuous interpolation between two physically distinct regimes: the deep ultraviolet (UV), 
where perturbation theory works well, and the deep infrared (IR), where non-perturbative 
techniques must be employed. In fact, the effective charge is intimately connected with two 
phenomena that are of central importance to QCD: asymptotic freedom in the UV, and 
dynamical gluon mass generation in the IR O]. Thus, while perturbatively it captures 
asymptotic freedom, it also exposes, due to the appearance of the Landau pole, the need of a 
non-perturbative regulating mechanism. Therefore, its low-energy behavior conveys valuable 
information about the way the theory cures the IR instabilities, namely through the non- 
perturbative generation of a dynamical mass scale, which tames the perturbative Landau 
pole. What makes the effective charge such an interesting quantity to study is its strong 
dependence on the detailed characteristics of some of the most fundamental Green's functions 
of QCD, such as the gluon and ghost propagators. Indeed, the basic ingredients that enter 
in its definition must contain the right information and be combined in a very precise way in 
order to endow the effective charge with the required physical and field-theoretic properties. 

In this article we will focus on two characteristic definitions of the effective charge, fre- 
quently employed in the literature. The first definition is obtained within the pinch technique 
(PT) framework , and its correspondence [7|, l8|] with the background-field method 

(BFM) jol. The PT effective charge, to be denoted by ap^{q'^), constitutes the most direct 
non-abelian generalization of the familiar concept of the QED effective charge. The second 
definition of the QCD effective charge, to be denoted by «gh(g^), involves the ghost and gluon 
self-energies, in the Landau gauge, and in the kinematic configuration where the well-known 



Taylor non-renormalization theorem 



Id, 



11| becomes applicable, agh(g^) has been employed 



extensively in lattice studies (see for instance [12, ll3l] and references therein), where the Lan- 
dau gauge is the standard choice for the simulation of the gluon and ghost propagators, as 
well as in various investigations based on Schwinger- Dyson equations (SDEs) IJ, ll5|]. Even 
though the theoretical origin of the two aforementioned effective charges is rather distinct, 
it turns out that, quite remarkably, by virtue of a powerful non-perturbative identity, they 
are almost identical in the entire range of physical (euclidean) momenta. In fact, most in- 
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terestingly, they are exactly equal in the deep IR {i.e., at vanishing momentum transfer, 

As we will see shortly, in the definition of the two effective charges appears a common 
ingredient, namely the gluon propagator (in the Landau gauge), and two ingredients that 
are not common. These two non-common ingredients are, a-priori, not related to each other; 
the role of the aforementioned identity is to furnish a non-trivial connection between them. 
Specifically, it relates the dressing function of the ghost propagator, denoted by F{q^), 
entering into the definition of aghl*?^), with a certain function, denoted by G{q^), appearing 
in the definition of a'PT{q^). The function G{q^) is the form-factor multiplying g^y in the 
Lorentz decomposition of a special Green's function, denoted by A^^{q), which appears 
in a variety of field-theoretic contexts. Most notably, A^^{q) enters in all "background- 
quantum" identities, i.e. the infinite tower of non-trivial relations connecting the BFM 
Green's functions to the conventional ones ly, ll7|. Notice also that G{q^) plays a central 



role in the new SDEs derived within the PT framework [18]; due to the special properties 
of the Green's functions involved, these new SDEs can be truncated in a manifestly gauge 
invariant way [sl. The identity in question connects the two non-common ingredients of the 
two charges, F{q^) and G{q^), to the second form factor of K^y{q), denoted by L{q^), in the 
way shown in Eq. (13.41) . 

To the best of our knowledge, the identity of Eq. (13. 4p was first derived in [lol, in con- 
nection with the so-called Kugo-Ojima confinement criterion [2^. The same identity was 



proved in 



211], where the general algebraic properties of SU{N) Yang Mills theories in the 



background Landau gauge were studied; however, no connection with the conventional 
Landau gauge was established. More recently, it was revisited in [22], where a new relation 
between the Kugo-Ojima parameter and the Gribov-Zwanziger horizon function has been 
advocated. However, to date, the dynamical equations for the quantities appearing in this 
identity remain largely unknown. 

In the present work we derive the central identity starting from the SDEs that govern the 
dynamics of the relevant functions, namely F{q^), G{q^), and L{q^). These SDEs allow for 
a detailed study of the individual properties of these three functions, both perturbatively 
and non-perturbatively. Most importantly, they expose the way these functions depend on 
the gluon propagator, and furnish a self-consistent framework for studying how an IR finite 
gluon propagator affects their IR properties. These properties, in turn, are responsible for 
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the mild discrepancy between the two effective charges mentioned above. 

The paper is organized as follows. In Section [Tll after introducing the necessary notation 
and definitions, we outline the basic theoretical ingredients entering into the construction of 
the two (dimensionful) renormalization-group (RG) invariant quantities, from which the two 
(dimensionless) effective charges, ap^i^q"^) and aghi^q"^) will be extracted. Then, we focus on 
the timely question of how to identify the correct non-perturbative scale in the presence of an 
IR finite gluon propagator. The central identity of the paper is derived in Section [ITI| starting 
from the defining SDEs. The renormalization procedure that preserves the validity of the 
identity is carried out, and various properties are studied in the UV and IR kinematic limits; 
most notably, we establish that if the gluon propagator is IR finite, then L{0) = 0. The 
implications of the identity on the two effective charges are discussed, and a relation between 
them is established, which is valid for the entire range of euclidean momenta. A detailed 
numerical analysis and comparison of the two effective charges at different renormalization 
scales is carried out in Section HVl using as an input the non-perturbative solutions of 
the SDEs corresponding to the various functions appearing in their definition. Finally, in 
Section |V] we present our conclusions. 

II. TWO NON-PERTURBATIVE EFFECTIVE CHARGES 

In this section we will first introduce some of the basic filed-theoretic ingredients necessary 
for the definition of the two effective charges we want to study. Then, we will briefly 
outline the basic construction and the assumptions involved in the definition of either charge. 
Finally, we will discuss in detail the important issue of how to extend the two definitions to 
the non-perturbative regime, and, in particular, the identification of the correct scale in the 
presence of an IR-finite gluon propagator. 

A. Definitions and ingredients 

Let us first introduce the notation and define some of the basic quantities entering into 
the problem under study. 

In the covariant renormalizable (R^) gauges, the gluon propagator A^^{q) has the form 
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FIG. 1: Diagrammatic representation of the functions H and A. 

where ^ denotes the gauge-fixing parameter, and P^v{q) = g^y — %(lv/<f is the usual trans- 
verse projector. Finally, A^^(g^) = + in(g^), with Ii^y{q) = P^j,(g)n(g^) the gluon 
self-energy (notice that since n(g^) has been defined with the imaginary factor i pulled out 
in front, it is simply given by the corresponding Feynman diagrams in Minkowski space). In 
addition, the full ghost propagator D{q^) and its dressing function F(g^) are related by 

D{q') = (2.2) 

Moreover, the all-order ghost vertex (after factoring out the color structure and the coupling 
constant g) will be denoted by r^(/c, q) with k representing the momentum of the gluon and 
q the one of the anti-ghost. The tensorial structure is given by 

- r^{k, q) = B,{k, q)q^ + B2{k, q)k^. (2.3) 

Thus, at tree-level r|f''(A;, g) = T^{k,q) = —q^j,. 

An important ingredient for what follows is the two-point function Af^^{q) represented in 
Fig. [H, defined by 

= -w'Ca [ Hf^^D{k + q)^P"{k)H,,{k,q), 
Jk 

= g,Mq')^'^-^L{q% (2.4) 

where Ca the Casimir eigenvalue of the adjoint representation [Ca = N for SU{Ny\, and 
= /i^'^(27r)~'^ j d'^k, with d = A — e the dimension of space-time. The scalar function G{q'^) 
appearing in the equation above allows the connection between the conventional and BFM- 
PT gluon propagators, and is known to play a central role in the PT formulation of the 
SDK. 
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The function Hf^^{k,q) (see Fig. [T] for a diagrammatic definition) is in fact a familiar 



object 



it appears in the all-order Slavnov- Taylor identity satisfied by the standard 



three-gluon vertex, and is related to the full gluon-ghost vertex by 

q''H^,{k,q) = -tT^{k,q). (2.5) 
At tree-level, Hju) = ig^j^u- Finally, using the most general Lorentz decomposition of if^^,, 

-iHf,y{k, q) = Ai{k, q)g^^ + A2{k, q)q^q^ + A^{k, q)k^ky + Ai{k, q)q^ky + A^{k, q)k^qu, (2.6) 
we obtain from (12.31) and (12. 5p two constrains for the various form-factors, namely 
Bi{k,q) = Ai{k,q) + q^A2{k,q) + {k ■ q)A4{k,q), 

B^ik, q) = {k- q)A,{k, q) + q^A^{k, q). (2.7) 



B. The pinch technique effective charge 

The QCD effective charges constructed within the PT uses QED as the basic refer- 
ence point [23]. In QED, one begins by considering the unrenormalized photon self- 
energy n°^(g) = g^P„/3(g)n°(g^), where Pap{q) = Qap ~ Qaqis/q^ and n°(g^) is a gauge- 
independent function to all orders in perturbation theory. After Dyson summation, we 
obtain the (process independent) dressed photon propagator between conserved external 
currents A|^^(g) = {gap/ q^)^^{q^)., with A°(g^) = + {q^)]"^ . The renormalization 
procedure introduces the standard relations between renormalized and unrenormalized pa- 
rameters: e = Z~^ef) = ZfZ]/"^ Z^^eo and 1 + in(g^) = Za[1 +in°(g^)], where Za (Zf) is the 
wave-function renormalization constants of the photon (fermion), Zi the vertex renormal- 
ization, and Ze is the charge renormalization constant. The Abelian gauge symmetry of the 
theory gives rise to the fundamental Ward identity (WI) 5'°T°(p,p+g) = Sq^ (p+q) — S^^ (p) , 
where r° and So{k) are the unrenormalized all orders photon-electron vertex and electron 
propagator, respectively. The requirement that the renormalized vertex r„ = ZiT^ and the 
renormalized self-energy S = ZJ^Sq satisfy the same identity, implies Zi = Zf, from which 

— 1/2 

immediately follows that Z^ = Z^ . Given these relations between the renormalization 
constants, and after pulling out the trivial factor gap/q^i we can form the renormalization 
group invariant combination, known as the effective charge, 

a(f) = ||A»tf ) = |iA(,^). (2.8) 
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In QCD, the crucial equality Z\ = Zf does not hold, because the Wis are replaced by 
;he more complicated Slavnov- Taylor identities (STIs), involving ghost Green's functions 
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23| . Furthermore, the gluon self-energy depends on the gauge-fixing parameter, already 
at one-loop order. These facts render the QCD generalization of a QED-like effective charge 
more complicated; however, the theoretical framework of the PT makes this definition pos- 



sible 



The PT rearranges the conventional gauge dependent n-point Green's functions, 
to construct individually gauge independent Green's functions, which, in addition, obey 
naive (ghost free) Wis . One important point, explained in detail in the literature, is the 

n n 

(all-order) correspondence between the PT and the Feynman gauge of the BFM [7|, |8| . In 



fact, using the methodology introduced in 2J], one can generalize the PT construction in 



such a way as to reach diagrammatically any value of the gauge fixing parameter of the 
BFM, and in particular the Landau gauge. In what follows we employ the aforementioned 
generalization of the PT, given that the identity we will eventually derive is valid only in 
the Landau gauge. 

The PT definition of the effective charge relies on the construction of an universal [i.e., 
process-independent) effective gluon propagator, which captures the running of the QCD f3 
function, exactly as happens with the vacuum polarization in the case of QED (See Fig. [2]). 
To fix the ideas, the PT one-loop gluon self-energy reads 



l + 6/ln(4 
/i 



(2.9) 



where b = llC^/487r^ is the first coefficient of the QCD /^-function. Due to the Abelian Wis 
satisfied by the PT effective Green's functions, the new propagator-like quantity A~^(g^) 
absorbs all the RG-logs, exactly as happens in QED with the photon self-energy. Then, the 
renormalization constants of the gauge- coupling and of the PT gluon self-energy, defined as 

A(g^/i2) = Z^\fi')Ao{q'), (2.10) 
where the "0" subscript indicates bare quantities, satisfy the QED-like relation 

Zg = Z-'/'. (2.11) 

Of course, Zg must be obtained under a given renormalization prescription, and the PT 
gluon self-energy will be then renormalized imposing (12.111) . Thus, regardless of the renor- 



(a) ib) (c) 

FIG. 2: The universal PT coupling, 
malization prescription chosen, the product 

^o(g') = gAiq') = g'if-^'Wq', i^') = d{q'), (2.12) 

retains the same form before and after renormahzation, i.e., it forms a RG-invariant (/i- 
independent) quantity 

For asymptotically large momenta one may extract from (i(g^) a dimensionless quantity 
by writing, 

where 'g'^{q^) is the RG-invariant effective charge of QCD; at one-loop 

^'^^'^ = TWM?M ^ Mn (gVA^en)' ^^'^^^ 
where Aqcd denotes an RG-invariant mass scale of a few hundred MeV. 

Eq. fl2.12l) is a non-perturbative relation; therefore it can serve unaltered as the starting 
point for extracting a non-perturbative effective charge, provided that one has information 
on the IR behavior of the PT-BFM gluon propagator A(g2). Interestingly enough, non- 
perturbative information on the conventional gluon propagator A(g2) may also be used, 
by virtue of a general relation connecting A(g2) and A(o^. Specifically, a formal all-order 



d[q') = (2.13) 



relation known as "background-quantum" identity 



16 



17| states that 



A{q')=[l + G{q')YA{q'). (2.15) 

Note that, due to its BRST origin, the above relation must be preserved after renormal- 
ization. Specifically, denoting by Zjy the (yet unspecified) renormahzation constant relating 
the bare and renormalized functions, Aq'^ and A^'^, through 

A'^'^(g,/.2) = ZA(/.2)Ar(g), (2.16) 



then from Eqs. (12.151) and (12. lip follows the additional relation 



(2.17) 



which is useful for the comparison with the coupling discussed in the following subsection. 

It is now easy to verify, at lowest order, that the 1+G{q^) obtained from Eq. (12.41) restores 
the {3 function coefficient in front of UV logarithm. In that limit [25] 



1 + 
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(2.18) 



Using Eq. (I2.15P we therefore recover the A~^(g2) of Eq. (12. 9p . as we should. 

Then, non-perturbatively, one substitutes into Eq. (I2.15P the 1 + G'(g^) and A(g2) ob- 
tained from either the lattice or SD analysis, to obtain A(g2). This latter quantity is the 
non-perturbative generalization of Eq. (12.90 : for the same reasons explained above, the com- 
bination 



d{q') 



[1 + G{q')] 



2 ' 



(2.19) 



is an RG-invariant quantity. 



C. The effective charge from the ghost-gluon vertex 

In the previous subsection it has become clear that the PT construction involves a par- 
ticular combination of two point functions only, with no explicit reference to any of the full 
vertices of the theory. Thus, as happens in QED, the effective charge so obtained is universal 
{i.e., it does not depend on the details of the process where the PT propagator is embedded), 
and depends naturally on a single scale, namely the physical momentum exchange of a given 
process. 

In principle, a definition for the QCD effective charge can be obtained starting from 
the various QCD vertices^, i.e., the ghost-gluon vertex, the three- and the four-gluon ver- 
tices, the quark-gluon vertex, etc ■ However, a priori, such a construction involves more 
than one scales, and further assumptions about their values need be introduced, in order 



^ In fact, as has been explained in detail in 26], an effective charge may also be defined from the gauge- 
invariant three-gluon vertex 
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to express the charge as a function of a single variable. As a general rule in all such a 
constructions one identifies a RG-invariant quantity formed by a judicious combination of 
the vertex form-factor and the self-energies associated with the fields entering into the ver- 
tex. Let us assume, for example, a vertex with three fields, i = 1,2,3, entering 
( with qi + q2 + Qs = 0). Denoting the corresponding propagators by Aj(gj), the relevant 
vertex form-factor by V{qi,q2,q3), by Zi the corresponding wave-function renormalization 
constants, and by Zy the vertex renormalization constant, one can renormalize the coupling 
such that^ Zg = Zy(ZiZ2Z^)~^^'^ , from which follows that the combination 

g2, gs) = 9^V\qi, q2, gs) Ai(gi) A2(g2) A3(g3) , (2.20) 

is a RG-invariant quantity. As mentioned above, the complication with this definition is 
that r^{qi,q2,q3) is a function of two kinematic variables. Thus, some additional assump- 
tion on the preferred kinematic configuration is usually introduced, such as, for example, 
?! = gi = gf = l'^ (and therefore qi ■ q2 = qi ■ qs = q2 ■ qs = — Q'^/2), which fully specifies the 
kinematic of the renormalization point. 

For the case of the ghost-gluon vertex, let us define in general the following renormaliza- 
tion constants 

A(g^/x2) = Z2\f^')A,{q'), 
F{q\^^^) = Z-\f^')Foiq% 
r^(A;,g,/x2) = Z^{fx'mk,q), 

go = Zg,{fi^)g'. (2.21) 

Notice that a priori Zgi defined as Zgi = ZiZj^^^'^Z^^, does not have to coincide with the Zg 
introduced in fl2.10p : however, as we will see in the next section, they do coincide by virtue 
of the basic identity we will derive there. 

In the Landau gauge, the form factor Bi of Eq. (12. 3p is UV finite at one-loop, and 
therefore, no infinite renormalization constant needs to be introduced at that order; of 
course, B2 must be UV finite in all gauges, and to all orders, otherwise the theory would be 
non-renormalizable. In order to obtain information about the UV behavior of Bi beyond 

^ In the MOM prescription, for instance, Zg is determined by requiring that the renormahzed vertex at the 
subtraction point assumes its tree-level value. 
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one-loop, one usually invokes the non-renormalization theorem of Taylor, which states that 
for vanishing ghost momentum (see Fig. [3]), one has that Bi{—q,q) + B2{—q,q) = 1, to 
all orders in perturbation theory. Given that B2 is finite to all orders (for any kinematic 
configuration), it follows that Bi{—q,q) is also finite to all orders. 

In particular, for the Taylor (vanishing incoming ghost momentum) kinematics, Zi will 
be determined as above explained by demanding that the relevant form factor be equal to 
its tree-level value after renormalization^, i.e., Z\ [(i?i(— g, g) + i?2(— g, g)] = 1. Then, one 
will have that 

Z\ = ZgiZj^ Z(^ = 1, (2.22) 

from which follows that 



Z-' = Z'fz,. (2.23) 



9 

Thus, the product 

f(g2) = /A(gV')i^'(gV') = goMq')F^iQ'), (2-24) 

forms either a dimensionful /x-independent combination or a UV cut-off independent one. 
Provided that we renormalize the propagators in the MOM scheme with Taylor kinematics 
(named as "Taylor scheme" in |l3]), /^(g^) is a RG-invariant combination. 

Therefore, for asymptotically large g^, in analogy to Eq. (12.131) one can define an alter- 
native QCD running coupling as 

f(g^) = (2.25) 
Notice that ^^/^(g^) has been shown to display the same behavior at any loop order as the 



ghost-gluon coupling for the Taylor kinematics (see Fig. [3]) in 13|. 
Using then Eq. (12.181) . and the fact that 



D'\q')=q 



9CAg\ fq^ 



IH ^ In , , 

4 487r 



(2.26) 



it is straightforward to verify that 'g^\i{q^) and 'g{q^) displays the same one-loop behavior, 
since, perturbatively the function 1+G'(g^) is the inverse of the ghost dressing function -F(g^). 
As we will see in the next section, this is nothing more than the one-loop manifestation of 
the more general identity relating G{q^) and F{q^). 



Recall that the form factor emerging at the Taylor kinematic limit fc^ — > — is Bi + B2. 
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FIG. 3: The ghost-gluon vertex and the Taylor kinematics. 
D. Effective charges from massive gluon propagators 

It is clear from the above analysis that before actually defining the effective charge with 
either method one constructs two dimensionful RG-invariant quantities, given in Eq. fl2.19p 
and Eq. (12.241) . with mass dimension -2. These two quantities share an important common 
ingredient, namely the scalar cofactor of the gluon propagator, A(g^), which actually sets 
the scale. The next step is to extract a dimensionless quantity, that would correspond to the 
non-perturbative effective charge. Perturbatively, i.e., for asymptotically large momenta, it 
is clear that the mass scale is saturated simply by g^, the bare gluon propagator, and the 
effective charge is defined by pulling a q^^ out of the corresponding RG-invariant quantity*^. 

Of course, as has been firmly established by now, in the IR the gluon propagator becomes 
effectively massive; therefore, particular care is needed in deciding exactly what combination 
of mass-scales ought to be pulled out. The correct procedure in such a case has been 
explained long time ago in the pioneering work of Cornwall and has been applied in 



various occasions 



]: a "massive" propagator, of the form [g^ + m^(g^)]^^ must be pulled 
out, where im?{q^) is a dynamical {i.e., momentum-dependent) mass^. 

Before applying this (correct) prescription to the two RG-invariant quantities in ques- 
tion, it is interesting to compare the situation with the more familiar, and conceptually 
more straightforward, case of the electroweak sector, where the corresponding gauge bosons 
{W and Z) are also massive, albeit it through an entirely different mass generation mech- 
anism. Specifically, while the W and Z bosons become massive at tree-level, through the 



^ This is equivalent to the standard MOM prescription for the coupling definition. 

^ Within the MOM philosophy one may implement the correct prescription by imposing 
A^^(/x^) = fj? + m?{(i'^) as the (non-perturbative) MOM renormalization condition for the gluon propa- 
gator. This prescription is equivalent to the standard one in the UV, while in the IR it introduces to the 
anomalous dimensions genuine non-perturbative (Borel non-analytical) terms of the type exp {—1/ gfi{q^)), 
which vanish as g'^ — > oo. 
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standard Higgs mechanism {i.e., fundamental scalars developing a vev), the gluons acquire 
their (momentum- dependent) masses non-perturbatively, through the dynamical realization 
of the well-known Schwinger mechanism {29I]. Despite the difference in their origin, the 
masses act in a very similar fashion at the level of the RG-invariant quantity associated 
with the corresponding gauge boson. 

Thus, in the case of the VT-boson, the corresponding quantity would read (Euclidean 
momenta) 



+ Ml 



with 



1 + &w5'w(/^) / dx\n 
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1 -1 



(2.27) 



(2.28) 



where bw = ll/247r^, and the ellipses denote the contributions of the fermion families. 
Clearly, rf^(O) = gl{0)/Ml, with ^^(0) = gi{fi)[l + b^gUfi) ln{Ml/ fi^)]~\ Evidently, in 
the deep IR, the coupling freezes at a constant value; Fermi's constant is in fact determined 
as 4:\^Gf = ^^^.(0)/M^. Note that in the case of QCD the corresponding combination, 
^^(0)/m^(0) would be similar to a Nambu-Jona-Lasinio type of coupling [30|]: at energies 
below the gluon mass m, the "tree-level" amplitude of four-quarks starts looking a lot like 
that of a four- Fermi interaction 3l[ |. 

This property of the "freezing" of the coupling can be reformulated in terms of what in 
the language of the effective field theories is referred to as "decoupling" 32|. At energies 
sufficiently inferior to their masses, the particles appearing in the loops (in this case the 
gauge bosons) seize to contribute to the "running" of the coupling. Possibly large logarithmic 
constants, e.g., ln(M|,//i^), may be reabsorbed in the renormalized value of the coupling. 
Of course, the "decouphng" as described above should not be misinterpreted to mean that 
the running coupling vanishes; instead, as already mentioned, it freezes at a constant, non- 
zero value. In other words: the "decoupling" does not imply that the theory becomes free 
(non-interacting) in the IR. 

This last clarification is not without relevance for the question at hand, namely the defi- 
nition of a physically meaningful effective charge. In particular, if one wants to extract an 



effective charge from an 
from lattice simulations 



R-finite g 



uon propagator (obtained from, e.g., SD studies [25| or 
35|), it would certainly be unwise to insist on the pertur- 
bative prescription, and simply factor out a l/g^. Even though one is merely redistributing 



33 



34, 
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RGI quantityl 




FIG. 4: The same RG-invariant quantity decomposed in two different ways, one giving a divergent 
propagator and a vanisliing coupling, and one giving a finite propagator and a finite coupling 

a given function, namely d^iq'^), into two pieces, factoring out deprives both of them 
of any physical meaning. Returning to the electroweak example, the effective coupling so 
defined would be given by the expression g^{q'^) = q^dw{q^), and so 5^(0) = 0; evidently, 
one would be attempting to describe weak interactions in terms of a massless, IR divergent 
gauge boson propagator and a vanishing effective coupling (See the curves in blue in Fig. HI). 
Given that the gluon propagator is finite in the IR, if this latter (wrong) procedure were 
to be applied to QCD, it would furnish a completely unphysical coupling, namely one that 
vanishes in the deep IR, where QCD is expected to be (and is) strongly coupled. 

As emphasized from the outset, the correct procedure is to factor out of the corresponding 
RG-invariant combination a "massive" propagator; in the FT case, we write the d{q'^) of 
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g2 _|_ rjj^2(^q2^ 

Given that d{q^) = g'^A{q'^), substituting Eq. fl^T^ into (^J^ we obtain 

c^..{g') = [<l' + m\g^)]^^^,, (2.30) 

where we have used apT^lq"^) = ^^(g^)/47r. As aheady mentioned, the dynamical mass m^(g^) 
appearing in the definition of a{q'^) is itself running; the explicit form of this running will 
be discussed in Section IIVI Similarly, from the RG-invariant quantity defined starting from 
the ghost-gluon vertex, given in Eq. (12.241) . we have that 

c^M) = «'(/^')(?' + m\q'))Aiq')F'{q'), (2.31) 

where a'{fi'^) = ^gh(/U^)/47r. 

Since A(0), F{0), G{0), and m(0) = rriQ are all finite (non- vanishing) , in the deep IR 
both couplings assume finite values given by 

agh(O) = mla'ifi^)A{0)[l+G{0)]-\ (2.32) 

III. DERIVATION OF THE IDENTITY FROM THE DYNAMICAL EQUATIONS 

In this section, we derive the central identity, valid only in the Landau gauge, relating the 
ghost dressing function with a particular combination of the form-factors G'(g^) and L[q'^) 
appearing in the tensorial decomposition of A^^, in Eq. (12. 4p . The proof hinges crucially 
on working in the Landau gauge (^ = 0), where the entire gluon propagator A^^{k) [and 
not just its self-energy n^j,(/c)] is transverse, i.e., k^^A^y{k) = 0. As we will see shortly, 
the operational consequence of this last property is that one can write q^A^y{k) = {q + 
k)^A^^{k), thus generating for free the appropriate ghost-gluon vertex, as needed. 

A. Deriving the relation 

The central relation is obtained as follows. First, consider the standard SD equation for 
the ghost propagator (Fig [5]), 

tD-\q^) =q^ + tg'CA [ T'^A^,{k)T''{k, q)D{q + k). (3.1) 
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k + q 



FIG. 5: The SDE for the ghost. 



Then, contract both sides of the defining equation fl2.4p by the combination g^g^ to get 



[G{q^) + L{q^)]q^ = g^CK [ q^^'"' {k) q'' H„,{k,q)D{k + q). 

Jk 



(3.2) 



Using Eq. (12. 5p and the transversaUty of the full gluon propagator, we can see that the rhs 
of Eq. (13. 2p is precisely the integral appearing in the ghost SDE (13. ip . Therefore 



[G{q') + W)]q' = iD-'\q')-q\ 
or, in terms of the ghost dressing function F(g^) [viz. Eq. (12.20 



l + G{q') + L{q') = F-\q'). 



(3.3) 



(3.4) 



The relation of Eq. ( 13. 4p . derived here from the SDEs of the theory, has been first obtained 
in 211], in the framework of the Batalin-Vilkovisky quantization formalism. As was shown 
there, the relation is a direct consequence of the fundamental BRST symmetry. 

Let us study the functions G{q^) and L{q^) more closely. From Eq. (12. 4p we have that 
(in d dimensions) 



1 



^ - > L{q') = JJZ:j^2 [dq^^'K^ ' ^'K) ' (3-5) 



{d- l}q 

which then gives, in terms of the SDE integrals 



d- 1 
d-l 



A'"'{k)H^p{k,q)D{k + q) + i\ [ q^A,^{k)T''{k,q)D{k + q) 

1 Jk 

4 / q'Ap,ik)T'^ik,q)Dik + q)+[ A^"^ {k) H,p{k,q)D{k + q) 
Q Jk Jk 



.(3.6) 



Inserting the decomposition of Eq. (12.30 and Eq. (12.60 into Eq. (13. 6p . and setting 

(k-qf 



fik,q)^ 



k'^q^ 



(3.7) 
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we obtain 

,2 



^^^2) ^)^^(^^ g) + [1 - /(A:, g)] g) - g^A^lA;, q)] } A{k)D{k + g), 

(3.8) 

while from Eq. (13.1 p 

F-\q') = 1 + s'Ca [[I- fik, q)]B,{k, q)^{k)D{k + g). (3.9) 
Jk 

Clearly, Eq. ( 13. 4p is automatically satisfied. 
B. Renormalization 

Of course, all quantities appearing in Eq. ( 13. Sp and Eq. ( 13. 9p are unrenormalized (we have 
suppressed the corresponding subscript "0" for simplicity); in particular, Eq. (13.40 involves 
unrenormalized G(g^), L^q"^), and F{q^). It is easy to recognize, for example, by substituting 
in the corresponding integrals tree-level expressions, that F~^[q^) and G{q^) have the same 
leading dependence on the UV cutoff Auv, namely 

i^u^(g^) = G,v(g^) = ^ln(^), (3.10) 

while Li^q^) is finite (independent of Auv) at leading order. The next step is therefore to 
carry out the necessary renormalization. 

As already mentioned above, the origin of the basic relation of Eq. ( 13. 4p is the BRST 
symmetry of the theory; in that sense, Eq. (13.40 has the same origin as the Slavnov- Taylor 
identities of the theory. Therefore, just as happens with the Slavnov- Taylor identities, 
Eq. (13. 4p should not be deformed after renormalization. Of course, the prototype example 
of such a situation are the Ward identities of QED; the requirement that the fundamental 
Ward identity g'^F^ = S^^ijp + g) — S^^ij)) should retain the same form before and after 
renormalization leads to the well-known textbook relation Zi = Z2 between the correspond- 



ing renormalization constants 23| . Similarly, for the case at hand, the renormalization must 
be carried out in such a way as to preserve the form Eq. (13.40 . Specifically, using the defi- 
nition given in Eq. (12.160 . in order to preserve the relation (13. 4p after renormalization, we 
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must impose that 

= Z,. (3.11) 

In addition, by virtue of (12. 5p . and for the same reason explained above, we have that, in 
the Landau gauge Ti,{k, q) and H„i,{k, q) must be renormahzed by the same renormahzation 
constant, namely Zi [viz. Eq. (I2.2ip ]: for the Taylor kinematics, we have that Zi = 1 [see 
Eq. i^M)]. 

Then, it is straightforward to renormalize Eq. (13. 8p and Eq. (13. 9p : using 



F-\q',f^') = Z,(A^„/i^)Fo-^(g^A^J, 
+ G(g^/i2) = Z,(A^^,^2)[l + Go(g^A^J], 

L(g2,/i2) = Z,{Al„f,')Lo{q\Al^), (3.12) 



we have that 



F-\q^) = Z, + g'CA [[I- f{k, q)]B,{k, q)/\{k)D{k + q) , (3.13) 

Jk 



and 



1+G(g^) = ^c+^ J{{d-l)A,{k,q) - [l-f{k,q)] [B,{k,q) - q'A2{k,q)]}A{k)D{k+q), 

(3.14) 

while the equation for Iv(g^) remains unchanged, i.e., one simply replaces in the second 
equation of (13. Sp the unrenormalized quantities by renormalized ones. This is consistent 
with the general observation made in 21|, according to which -L(g^) does need its own coun- 
terterm, i.e., one proportional to q^qu, in order to get renormalized. The situation is similar 
to what happens with the CT/^i^q" part of the standard QED vertex: The renormalizability 
of the theory forbids of course a counterterm proportional to such a tensorial structure; the 
magnetic form factor (usually denoted by F2{q^)) is made finite (beyond one loop) after 
multiplication by the renormahzation constant Zi (whose counterterms are proportional to 
7^). Thus, while the one-loop answer for F2 is finite, at higher orders one gets divergences 
proportional to cr^^q^ which are, however, canceled exactly (order by order) by the inclusion 
of the Zi counterterms in the Feynman graphs of the previous order. For this reason, just as 
F2, despite its one-loop finiteness L depends in general on the UV cutoff A^^, as indicated 
explicitly in Eq. (137721) . 
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C. Calculations and approximations 



In order to study the relevant equations further, we will approximate the form fac- 
tors Ai{k,q) and Bi{k,q) with their tree-level values, i.e., Ai{k,q) = Bi{k,q) = 1, and 
A2{k, q) = 0; according to lattice studies 3^, this appears to be a very good approximation. 
Then, we obtain from Eqs. (13.81) and (13.91) 

F'\q') = Z^ + g'C^ [[l~f{k,q)]A{k)Dik + q), 

Jk 

1 + G{q') = Z, + ^lj,{d-2) + f{k,q)]A{k)D{k + q), 

^^^') = f-IT [[^-dfik,q)]Aik)Dik + q). (3.15) 
Now, it turns out that if F and A are both IR finite, then 



[l-dfik,q)]Aik)D{k + q) 



= 0, (3.16) 



To see this, one may use the result k^kyF{k)A{k) = g^yd ^ F{k)A{k), or, equivalently, 
go to spherical coordinates and use that® 



d^sin'*^(l-dcos'^) = 0. (3.17) 
Thus, from Eq. (I3.15P we obtain the important result 

L(0) = 0, (3.18) 
under the assumption that F and A are IR finite. In addition, using ( 13.160 . we obtain 

F-i(O) = 1 + G'(O) = Z,+ ^'^a(c/-1) I ^(j^^jjfj^y (3_^g) 

" Jk 

Note that perturbatively, at one loop, Eq. (I3.16P does not hold, because in that case 
A{k) is not IR finite; consequently, at one loop L(0) ^ 0. Specifically in this case, using 
dimensional regularization, we obtain the g-independent result 

l~df{k,q) 3 I 



k k^{k + q)^ 2 167r2 



(3.20) 



6 Recall that Cde sin" 9 = ^^"fllffl^ 

H 2 j 
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which gives 



9^ 
327r2 ■ 



(3.21) 



If instead we were to use an IR finite gluon propagator, modeled simply by 
A^^(A;) = k"^ — m?, the same calculation would show that Lm{q^) depends non-trivially on 
[see Eq. (13.261) below], and in fact, -Lm(O) = 0. 

We next go to the Euchdean space, by setting —q^ = g|, and defining AE(g|) = — A(— g|), 
De(9e) — ~-^(~5'e); for the integration measure = i ff. ■ Then, using Eq. (12. 2p and 



suppressing the subscript "E" , we obtain from Eqs. (13.81) and (13.91) 

F~\q') = Z,-/Ca l[l-f{k,q)]^{k)D{k + q), 



1 + G{q^) 
L{q') 



9'Ca 
d~ 1 



[{d-2) + f{k,q)]A{k)D{k + q), 



[l-df{k,q)]A{k)Dik + q). 

, rv 



(3.22) 



d — ^ Jk 

Next let us introduce spherical coordinates. Setting = x, k'^ = y, we have that 
k ■ q = ^Jxy cos d, and so ik ■ qf jq^ = y cos^ 9, and {k + qy = x + y + 2y/xy cos 9. Moreover, 
at d = 4, the measure is given by 



d'^k = 27r I d9 sin^ 9 dy y. 
Jo 



(3.23) 



Let us first consider the case in which the ghost propagator assumes its tree-level form, 
namely D{k + q) = l/{k + g)^. Then, using the results 
sin^^ 



d9 



d9 



X + y + 2^/xy cos 9 
sin^ 9 cos^ 9 



2 



^0 X + y + 2.^/xy cos9 
where G(x) is the Heaviside step function, one obtains 



— 6(x — y) -\ — Q{y — x) 
X y 

^(l + y)Q(x-y) + -(l + -]Q{y-x) 
x\ xj y\ y' 



(3.24) 



1 + G{x) 
L{x) 
F-Hx) 



1271 



^ j\ y[^ + £) A(z/) + jjy (3 + ^) A(y) 

1 rdyy'A{y)+x Hdy^ 

Jo Jx y 



IGvr 



- fdy yh--) A{y) + Hdy (3 - -) A{y)\ . (3.25) 

X Jo V xJ V yJ 



Substituting into the equation for L{x) the tree-level value for A{y) we obtain the constant 

result Lo{x) of Eq. (I3.2ip . On the other hand, using A{y) = {y + m^)~^, we find 

OsCa f 1 Ix"^ 2 4 , / ^ 
[X] = —rz — <— — — mx + m In (IH 



12n 



x" 



X I m 
+ — In 1 + — 



X 



(3.26) 
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from which we clearly see that Lm{0) = 0. In addition, for large x, L„i{x) goes over to the 
massless limit of Eq. (13.211) . 

The general case for an arbitrary ghost dressing function F{k + g), can be treated by 
means of the angular approximation. Specifically, one write approximately 



1 + G{x) 
L{x) 
F-Hx) 



12n 



IQlT 

F{x) 



F{x) 

X 

X 



X^ 



f\ y (3 + ^) A(y) + jjy + ^ A{y)F{y) 

px poo 

I dy y'^A{y) + x dy 

Jx 



A{y)F{y) 



IGtt 



F{x) 



X 



/X /*oo / 

dyy{3~ ^) A{y) + j dy [?, ~ 



A{y)F{y) 



(3.27) 



It is then easy to see {e.g.^ by means of the change of variables y = zx) that if A and F 
are IR finite, then L(0) = 0, as claimed before. Let us now assume that the renormalization 
condition for F{x) was chosen to be -F(/i^) = 1. This condition, when inserted into the third 
equation of ( ]3.27p . allows one to express Zc as 



IGtt 



P Jo 



dyy 3 



A{y) 



3 - 



A{y)F{y) 



(3.28) 



and may be used to cast (13.271) into a manifestly renormalized form. Note that if one choses 
-F(/i^) = 1 then one cannot choose simultaneously G{^'^) = 0, because that would violate 
the identity of Eq. (13. 4p . given that L(/i^) 7^ 0. In fact, once F{fi^) = 1 has been imposed, 
the value of G{^^) is completely determined from its own equation, i.e., the first equation 
in ([32ZD. 

In addition in the MOM scheme the conventional and PT propagator can- 
not be made equal at the renormalization point, since the identity (12.151) implies 
AM)(;.^)=^2[l + G2(/.^)]l 



D. Implications for the effective charges 

After this general discussion, let us now return to the couplings, and discuss the impli- 
cations of the identity and the dynamics we have derived. 

First of all, comparing Eq. (I2.12p and Eq. (I2.24p . it is clear that g{jji) = g'{fi), by virtue of 
Eq. (13.111) . Therefore, using Eq. (I2.15p . one can get a relation between the two RG- invariant 
quantities, r(g^) and d{q'^), namely 

r{Q') = [^ + G{q')rF\q')d{q'). (3.29) 
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From this last equality follows that Opx and agt{q'^) are related by 



a,^{q') = [l + Giq')YF\q')aMq 



(3.30) 



After using Eq. (13.41) . we have that 

«gh(g^) = 

or, equivalently, 



1 + G(g2 



aprr{q ) = agh(g ) 



1 + G(g2) + L(g2) 

Hq' 



apT^q^). 



1 + 



1 + G(g2) 

Evidently, the two couplings can only coincide at two points: (i) at 
the fact that L(0) = [see Eq. (I3.18p ]. we have that 

agh(O) = apT(O), 



(3.31) 
(3.32) 

0, where, due to 
(3.33) 



and (ii) at = oo, given that in the deep UV L^q"^) approaches a constant. Note in fact 
that the two effective charges cannot coincide at the renormalization point n, where 



agh(/i ) = [1 - LifJ' )] apT(/i ) 



(3.34) 



this can be understood also in terms of the discussion following Eq. f l3.28p . 

As we will see in the next section, the numerical analysis reveals that L{q'^) is fairly small 
compared to G{q^); thus, even in the region of intermediate momenta, where the difference 
reaches its maximum, the relative difference between the two charges is less than 5%. 



IV. NUMERICAL ANALYSIS 



In this section we will compute the QCD effective charges defined above, using as input for 
the various Green's functions appearing in their definitions the non-perturbative solutions 
of the corresponding SDEs, in the Landau gauge. In particular, we will solve numerically a 
system of three coupled non-linear integral equations, containing A(g^), -F(g^), and G{q^) 
as unknown quantities. Once solutions for these three functions have been obtained, then 
L(g^) is fully determined by its corresponding equation, namely the second one in Eq. (13.271) . 
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FIG. 6: The new SDE for the gluon propagator at the one-loop dressed level. 



A. The system of SD equations 

The two SDEs determining -F(q^) and G{q^) are given in Eq. fl3.27p . The SD equation 
governing A(g^), is given by 25| 



[1 + G{q^)f/\-\q^)p,M = q'PM + ^ E' 



(4.1) 



1=1 



where the diagrams {ai)^i, are shown in Fig. [61 As explained in 25|, due to the abelian 
Ward-identities satisfied by the fully-dressed vertices in the PT-BFM scheme, we have that 
(l^[{C'i)fMu + {(^2)fiu] = (l'^'[{(^3)i^iu + ('^4)^4!^] = 0. This last property enforces the transversality 
of the gluon self-energy "order-by-order" in the dressed-loop expansion, which is one of the 
central features of the gauge-invariant Schwinger- Dyson truncation scheme defined within 



the PT-BFM framework 



18| 



After introducing appropriate Ansatze for the aforementioned fully-dressed vertices, we 
finally arrive at the integral equation 



+ 9'Ca 



9'Ca 
6 

-/ 



Jk Jk ^ Jk 



k^{k + qY 



D{k)D{k + q)-2 / D{k) 

'k 



(4.2) 



with 



/i = 20g2 + 18A;2 - 6(A; + g)2 + 
, k^ 



2\2 



{k . qf 



_ 27 

/2 - -y 



+ 



{k + qf 
2 {k-qf 



+ 4- 



20 10 g2 
\ \ 

j/2 g2 k'^[k + qY 

{k-qf 



2{k + qy 
q2^2 



-4 



(4.3) 



{k + qf {k + qy k^{k + qY q^{k + q^ 

The important point is that, by virtue of the massless composite poles introduced into the 
SDE through the particular Ansatze employed , one obtains an IR finite solution 

for the gluon propagator, i.e., a solution with A^^(O) > 0, in complete agreement with a 



large body of lattice data 
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35|. As explained in detail in |25|, the formal expression 
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Ghost dressing 

a(M') = 0.21 and n = 4.3 GeV 
— ^a{^i') = 0.16andn= 10 GeV 
— afr') = 0.13andn = 22 GeV 



1E-4 1E-3 0,01 0,1 



100 1000 




1E-4 1E-3 0,01 



q'[GeV'] q'[GeV'] 

FIG. 7: Left panel: Numerical solutions for the gluon propagator obtained from the SDE using 
three different renormalization points: fj, = 4.3 GeV and a{fj,'^) = 0.21 (black curve), fj, = 10 GeV 
and a^fi^) = 0.16 (red curve), /i = 22 GeV and a(;U^) = 0.13 (green curve). Right panel: The ghost 
dressing function F{q'^) obtained from its corresponding SDE and renormalized at the same points. 



determining A^^(O) involves quadratically divergent integrals, which may be regulated using 
the standard rules of dimensional regularization. This procedure leaves the (finite) value of 
A^^(O) largely undetermined; therefore, in practice, A^^(O) is treated as a free parameter, 
whose value is to be fixed using phenomenological constraints or lattice data. In addition, 
and because A~^(0) is finite, the ghost dressing function F[q^) clearly saturates in the deep 
IR, reaching a finite value at = (no "enhancement" observed), in agreement with recent 



lattice data 



35l |. and a variety of independent studies [15 



B. Solutions and checks 

In Fig. [71 we show the numerical results for A(g^) and F{q'^), renormalized at three 
different points. On the left panel, the black curve represents the numerical solution of 
A(g^) when a(/i^) = 0.21 and fi = 4.3 GeV. The red curve is obtained when a(/i^) = 0.16 
and /i = 10 GeV, while for the green curve we used a(/i^) = 0.13 and /i = 22 GeV. On the 
right panel we plot the corresponding F{q^) renormalized at the same points. 

In Fig. [8] we show the numerical results for the functions 1 + G{q^) and L{q'^), using 
the same renormalization points used previously. The color pattern is also the same as 
before. For values of < O.lGeV^, we then see that [1 + G'(g^)]^ develops a plateau and 
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q'[GeV'] q'[GeV'] 



FIG. 8: Left panel: 1 + G{q^) determined from Eq. ([3221), using the solutions for A(g^) and D{(f') 
presented in the Fig.[7]at the same renormahzation point. Right panel: The function L{q^) obtained 
from Eq. (13:271) . 

saturates at a finite value in the deep IR region. In the UV region, we instead recover the 
perturbative behavior f l2.18p . On the other hand, L{q'^) (right panel) shows a maximum in 
the intermediate momentum region, while, as expected, L{{)) = 0. 

With all ingredients defined, the first thing one can check is whether Eq. 02.191) gives rise 
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FIG. 9: The product d{q'^) obtained combining the results for A(g^) and [1 + according to 

Eq. dllSI). 
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FIG. 10: The product (1 + G)F for different values of the renormalization point fi. Note the fine 
scale of the y-axis. 

to a RG-invariant combination, as expected. Using the latter definition, we can combine the 
different data sets for A(g^) and [1 + G'(q'^)]^ at different renormaUzation points, to arrive 
at the curves shown in Fig. [91 Indeed, we see that the combination d{q'^) is practically 
independent of the renormalization point chosen. 

In addition, from the available solutions we can compute the product {1 + G)F, which, ac- 
cording to Eq. (13.301) . relates the two effective charges of interest. Evidently, since both effec- 
tive charges are supposed to be RG-invariant quantities, so should be the product (1 + G)F 
relating them. In Fig. [10] we plot {1 + G)F for different values of the renormalization point /x; 
clearly the dependence on /x is very mild. The theoretical origin of this residual /i-dependence 
can be traced back to the approximations used for the ghost-gluon vertex T,^ and the func- 
tion if^jy (see beginning of subsection C). This approximation distorts the multiplicative 
renormalizability of the corresponding SDEs; indeed, for multiplicative renormalizability to 
be enforced, one must assume the exact renormalization properties for and H^,^, as was 
done in subsection B, where the renormalization was carried out formally. Instead, the ap- 
proximation employed causes a mismatch in higher orders, which introduces the observed 
mild dependence on fi. This dependence can be eliminated by resorting to the systematic 
improvement of the corresponding Ansatz used for T,y, in the spirit of the prototype QED 
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FIG. 11: The behavior of the running mass given by Eq. ()4.4I) when mo = 500 MeV (blue hue) and 
uiQ = 600 MeV (magenta line). In both cases we used Aqcd = 300 MeV. 



calculations presented in 39|], and more recently in 40 1. 



C. The effective charges 

We can next proceed to extract the non-perturbative running charge ap^lq"^), defined in 
Eq. fl2.30p . by multiplying the results obtained for d{q'^) by the factor [g^ + ^^( q^)] - To this 
end, we will assume that 7ii?{q^) has a power-law type of running, given by 

I 3 
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(4.4) 



A variety of theoretical and phenomeno 
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43j. In Fig. Owe 



q'^ + mg 

Notice that when ^ one has m^(0) 
logical estimates place it in the range m^ = 350 — 700 MeV 
plot the behavior of m^{q^) as given by Eq. (14. 4p . for the two values mo = 500 MeV and 
mo = 600 MeV, which will be used in the rest of this section. 

On the left panel of Fig. [121 we show the results for apT(g^) when mo = 500 MeV in 
Eq. (14. 4p . The small discrepancy between the three curves is mainly due to the propagation 
of the tiny residual /x dependence displayed by the quantity d{q^) as shown in Fig. [9l One 
clearly sees that the effective coupling a^^^q^) freezes out and acquires a finite value in the 
IR, while in the UV it shows the expected perturbative behavior. For mo = 500 MeV, one 
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FIG. 12: Left panel: The running charge obtained from (j2.30p using the SDE solutions for A(g^), 
D{(f'), and \ + G{(f'). We use a running mass given by Eq. ()4.4p with ttiq = 500 MeV. Right panel: 
The same for rriQ = 600 MeV. 

gets apT(O) ~ 0.6. One should also notice that the choice of smaller values of mo would 
not produce a monotonically decreasing apT(g^); instead, one observes the appearance of 
"bumps" in the IR region. Therefore if one were to introduce the monotonic decrease as 
an additional requirement of ap^i^q^), this would provide a lower bound for the possible 
values of mo. Finally, on the right panel of Fig. [121 we show the effective coupling for the 
case mo = 600 MeV. Now, the freezing occurs at the slightly higher value of apT(O) ~ 0.85. 
Evidently, the freezing value apT(O) increases as one goes to higher values of mo. 

An accurate fit for the running charges shown in Fig. [12] is provided by the following 
functional form 
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with the function h{q ,m (q )) given by 



h{q\ m\q')) = p,m\q^) + p.^^^^y (4.6) 

Our best fits to the numerical results for ap^lq"^) using Eq. (14. 5 p above are shown in Fig. [131 
Finally, we compare numerically the two effective charg (g^) and agh('?^)- The 

results are shown in Fig. [HI where r(g^) is compared with d{q'^) (left panel), and agh((3'^) 
with apT(g^) (right panel). As anticipated, the curves coincide in the deep IR and UV, and 
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FIG. 13: Left panel: The fit given by Eq. (|4.5p for niQ = 500 MeV; in this case the best fit values 
correspond to pi = 4.5, and p2 = —2. Right panel: Same as before in the case mo = 600 MeV; in 
this case the best fit parameters are pi = 2.2, and p2 = —1.25. 

differ only slightly in the intermediate region. To produce both curves, we have factored out 
a mass of mo = 500 MeV, whose dynamical running is again given in Eq. (14.41) : equivalently, 
one could use directly Eq. f l3.30p . 



V. CONCLUSIONS 



In this article we have presented a detailed comparison between the two QCD effective 
charges, apT(g^) and aghil"^), obtained within two vastly different frameworks: the PT (and 
BFM) on the one hand, and the ghost-gluon vertex (with the Taylor-kinematics) on the 
other. It turns out that their dynamics involves the gluon propagator A(g^) (in the Landau 
gauge) as a common ingredient, entering in both apT^q"^) and aghl?^), and two different 
ingredients, which participate in a non-trivial identity. This identity, which is valid only 
in the Landau gauge, relates the ghost dressing function, -F(g^), with the two form-factors, 
G(g^) and L{q^), appearing in the Lorentz decomposition of a special Green's function, 
originating from the ghost sector of the theory. 

The two QCD effective charges have been computed using as input the non-perturbative 
solutions of a system of three coupled non-linear integral equations, first derived in 25| . 
containing A(g^), F{q^), and G{q^) as unknown quantities. The solutions obtained from the 
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q'[GeV'] 



q[GeV] 



FIG. 14: Left panel: Comparison between the RG- invariant products r(g^) (red curve) and d{q'^) 
(black curve). Right panel: Comparison between Qgh((7^) (red curve) and ap^iq^) (black curve) 
when a mass of niQ = 500 MeV is factored out. In both plots the difference between the curves 
appear in the intermediate regime of momenta being entirely due to the function L(g^). 

above system of SDEs for A(g^) and -F(g^) - and subsequently fed into the defining equations 
of the effective charges- are in qualitative agreement with recent results from large- volume 



lattices, both for SU{2) [3J] and SU{3) [35|: both quantities reach finite (non- vanishing) 



values in the deep JR. One important consequence of the central identity (and the dynamics 
encoded in the relevant equations) is that the two charges are identical not only in the deep 
UV, where asymptotic freedom manifests itself, but also in the deep IR, where they "freeze" 
at the same non- vanishing value. 

As already mentioned in section \IV\ at the level of the SDE for the gluon propagator, 
namely Eq. 04.21) . the value of A(0) is a free parameter. The value chosen for A(0) affects (in 
a non-linear way) the IR values of the RG-invariant quantities, namely d{0) and r(0), which, 
in turn, restricts the values of the gluon mass, mo, and the freezing value of the effective 
charges. Throughout the analysis presented in section [IV] the criterion used for choosing 
the values of A(0) was that the resulting values for mo and agh(O), (or apT(O)) would be 
numerically compatible with those obtained from a variety of phenomenological studies 43|. 
Specifically, values for mo in the range of 350 — 700 MeV and q;pt(0) ~ 0.7 ± 0.3. Notice, 
however, a subtle point that may be of relevance when carrying out such comparisons. 



30 



The effective charge assumed in most of the aforementioned studies is that of [ij, which 
has a very particular functional form, and corresponds to the standard PT construction, 
where the Feynman gauge of the BFM is dynamically singled out. Instead, for the reason 
explained in subsection 2.2, the present analysis is based on the generalized PT [24], which 
eventually projects one to the Landau gauge of the BFM. It would be interesting to reach 
a quantitative understanding of how the aforementioned difference in the gauges affects 
the phenomenological values of the gluon mass and of the freezing of the effective charge. 
Calculations in this direction are already in progress. 



As has been emphasized in [25|], even though the solutions of the SDE system are in 
qualitative agreement with the aforementioned lattice results, they display a considerable 
quantitative discrepancy from them. Specifically, A(g^) differs significantly in the region of 
intermediate momenta, and the value of the ghost dressing function is about a factor of two 
less than that obtained on the lattice. These discrepancies, in turn, are expected to affect 
the numerical values (but not the qualitative features) of quantities computed using them 
as input. In particular, it should be interesting to obtain the QCD effective charges studied 
here using as input the lattice results for A(g^) and F{q^), and [indirectly, using, e.g., the 
first equation in fl3.27p ] for G{q^)\ we hope to address this issue in a future work. 
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